Quasi-conservation laws for compressible 3D Navier-Stokes flow 
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We formulate the quasi-Lagrangian fluid transport dynamics of mass density p and the projection 
q = u • Vp of the vorticity oj onto the density gradient, as determined by the 3D compressible 
Navier-Stokes equations for an ideal gas, although the results apply for an arbitrary equation of 
state. It turns out that the quasi-Lagrangian transport of q cannot cross a level set of p. That is, in 
this formulation, level sets of p (isopychnals) are impermeable to the transport of the projection q. 
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PACS numbers: 47.10.A-, 47.10.-g 

The aim of this note is to formulate the fluid conser- 
vation dynamics for mass density p and the projection 
q = ui ■ Vp of the vorticity uj = curlw onto the mass 
density gradient, as determined by the 3D compressible 
Navier-Stokes equations [1] 

Du 

p-^— =/xAit — Vo7, w = p — (p/3 + p")div u , (1) 



— +pdivw = 0, — =d t +u-y, 
DO p 

c v — = - div u + Q , p = RpO . 
Dt p 



(2) 
(3) 



Here, it denotes the spatial fluid velocity, p is the shear 
viscosity and p v is the volume viscosity, both of which 
are taken as constitutive constants of the fluid. 

For definitcness, we have chosen an ideal gas equation 
of state to relate pressure, p, temperature, 6, and mass 
density, p. In addition, R is the gas constant, c v the 
specific heat is constant and Q is the heating rate, which 
we may assume is known. The pressure depends on the 
two thermodynamic variables p and 8. It is noteworthy 
that the dynamics of the temperature and the choice of 
equation of state in (3) will not affect our considerations 
below of the transport dynamics of the projection q = u>- 
Vp. This means the geometric considerations that follow 
are universal for any viscous compressible fluid flow. 

According to equations (1) and (2) the vorticity u> = 
curl u evolves according to a stretching and folding equa- 
tion 



dt<jj — curl (u x u>) 

= np~ 1 Auj + Vp -1 x 



(4) 



pAu — V ( TX7 — — 



driven by the right hand side. The form of the last term 
invites a projection against Vp. By the product rule, 
the projection q = u> ■ Vp satisfies a geometric relation 
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reminiscent of Ertcl's Theorem [2] 

Dq /Du „ \ „ „ [Dp\ 

— - = w • Vm • Vp + u> • V — 

Dt \Dt / V Dt J 



(5) 



Ertel's Theorem, important in atmospheric dynamics, 
has the same form as (5), with mass density replaced 
by potential temperature. Substituting mass conserva- 
tion (2) and Navier-Stokes vorticity dynamics (4) into 
the purely geometric equation (5) yields the following 
equation for q written in a suggestive divergence form, 

d t q + div qu + div(u;pdivM — pAu x V(ln p)j = . (6) 

Now we apply an observation of Haynes and Mclntyrc 
[4] that first arose in atmospheric physics and allows one 
to define a current density J. Assuming that solutions 
exist for equations (l)-(3), as 

J = qu + up div u — pAu x V(ln p) + V</> x V/(p) , (7) 

where <j> is an undetermined gauge potential and / is an 
arbitrary diffcrentiable function of p. Then (6) and (2) 
can be rewritten in the quasi- conservative form 



dtq + div J = and qdtp + J • Vp = . 



(8) 



The relation J • Vp = q div pu allows zero projection, 
q = 0, by the second equation in (8). Thus, the projection 
q may vanish anywhere in the flow, but it cannot be 
maintained, because div J ^ 0. Together, the equations 
in (8) imply a family of conserved quantities, since 



ft(g*'(p))+div (J4>'(p)) = 0. 



(9) 



for any function $'(p) = d^/dp. 

The conserved densities q$'(p) = div($(p)u;) in (9) 
possess quite different flow properties from those of mass, 
energy and momentum. The question of the physical 
interpretation of these quantities may now be raised. 
We know that equations (8) are purely kinematic, be- 
cause the projection taken against Vp removed any de- 
pendence on the dynamics of the temperature and the 
choice of equation of state in (3). Nonetheless, equa- 
tions (8) have been derived without approximation from 
the Navier-Stokes fluid equations for compressible motion 
and mass transport. 
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Moreover, their analogues also occur and have been 
found useful in other areas of fluid dynamics, particularly 
in the dynamics of the ocean and atmosphere [3, 4]. (In 
the oceanic context, p denotes buoyancy and the motion 
is usually taken to be incompressible.) 

Indeed, a similar calculation may be performed to de- 
rive the quasi-conservative form in (8) of the dynamics 
of the projection of vorticity on temperature gradient, 
q' = ui ■ V8 from equation (5) with V p replaced by V0. 
The projection q' with 6 denoting the potential temper- 
ature is known as the potential vorticity density and has 
been useful in guiding thinking about nonlinear convec- 
tive processes in atmospheric science, [3, 4]. 

In addition, the corresponding forms of (7)-(8) may be 
written in the absence of viscosity. One hopes that with 
the proper interpretation equations (8) may also be useful 
in considerations of compressible Navier- Stokes flows. 

Following [3, 4] in the atmospheric context for which 
the projection q is replaced by potential vorticity q' , we 
propose to interpret equations (7)-(8) locally as a type of 
impermeability theorem for the quasi- Lagrangian trans- 
port of the projection q and the mass density p by a 
pscudo- velocity U = J /q. Here, impermeability means 
that the projection q = u) ■ Vp cannot be transported by 
pseudo- velocity U [7] across level sets of the mass density, 
p. This is remarkable because the two quantities q and p 
are transported by different velocities. 

The quasi-conservative form of compressible Navier- 
Stokcs fluid dynamics in (8) has apparently not appeared 
previously in the literature, and we hope that this formu- 
lation may find some future use, perhaps in analogy with 
the use of potential vorticity in the atmospheric context. 

In particular, these equations may be useful in the 
study of stretching and folding in compressible fluid 
flows, just as has been recently investigated in the at- 
mospheric context in [5, 6]. Although previous studies 
of compressible flows have often focused on shock forma- 
tion, we shall close this note with a remark about how the 
gradient Vq of the projection q = u) ■ Vp participates in 



stretching, folding and expansion of higher-order gradi- 
ents in compressible Navier-Stokes fluid flows. In prepa- 
ration, we rewrite Equations (2) and (8) in terms of the 
pseudo- velocity U = J/q as 

d t q + div (qU) = d t p + U • Vp = . (10) 

Note however that divW ^ 0. We define 

S=VgxVp, (11) 

and discover from a direct computation that B satisfies 

d t B-cm\(UxB) =T>, (12) 

where T> = —\7(qdivU) x Vp. The proof of the corre- 
sponding relation in an atmospheric physics context can 
be found in [5, 6]. 

Based on the transport velocity defined in (7), the left 
hand side of (12) makes it clear that the vector B under- 
goes the same type of stretching, folding and expansion 
processes driven by the X>-vector on the right hand side 
as occurs in the vorticity equation (4). Explicitly, the 
stretching, folding and expansion processes are given by 

d t B+ U-VB-B-VU + B&ivli = T> . (13) 
" * ' 

stretch, fold & expand B 

This is remarkable, because B = Vq x Vp contains infor- 
mation not only about Vu> but also about Vp and even 
VVp. Thus, the stretching and folding in the original 
vorticity equation (4) compounds itself in the same form 
in the S-equation (12), but with higher spatial deriva- 
tives. 
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